Quantifying modeling uncertainties when combining multiple gravitational-wave detections from binary neutron star sources

With the increasing sensitivity of gravitational-wave detectors, we expect to observe multiple binary neutron-star systems through gravitational waves in near future. The combined analysis these signals offers possibility constrain radius and equation state dense nuclear matter with unprecedented accuracy. However, it is crucial ensure that uncertainties inherent models will not lead systematic biases when information from detections combined. To quantify waveform systematics, perform an extensive simulation campaign sources analyze them a set four different models. For our 38 simulations, find statistical decrease $\ifmmode\pm\else\textpm\fi{}250\text{ }\text{ }\mathrm{m}$ (2% at 90% credible interval) but differences between currently employed can be twice as large. Hence, essential become dominant inferences capitalizing on future developments.